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ABSTRACT 

The development of fully nonlinear Gortler vortices in the high Reynolds number flow 
in a symmetrically constricted channel is investigated. Attention is restricted to the case of 
‘strongly’ constricted channels considered by Smith and Daniels (1981) for which the scaled 
constriction height is asymptotically large. Such flows are known to develop a Goldstein 
singularity and subsequently become separated at some downstream station past the point 
of maximum channel constriction. It is shown that these flows can support fully nonlinear 
Gortler vortices, of the form elucidated by Hall and Lakin (1988), for constrictions which 
have an appreciable region of local concave curvature upstream of the position at which 
separation occurs. The effect on the onset of separation due to the nonlinear Gortler modes 
will be discussed. A brief discussion of other possible nonlinear states which may also have 
a dramatic effect in delaying (or promoting) separation will be given. 


1 Research was supported by the National Aeronautics and Space Administration under NASA Contract 
No. NASM8605 while the authors were in residence at the Institute for Computer Applications in Science 
and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1. Introduction 


Our concern is with the development of fully nonlinear Gortler vortices in Smith- 
Daniels flows and their effect on the onset of separation. By Smith-Daniels flows we 
refer to the steady flow in a symmetrically constricted channel considered by Smith and 
Daniels (1981) who showed that when the scaled height, h 1, of the constriction becomes 
large a secondary classical boundary layer of thickness 0{h ~ a / 2 ) is set up within the wall 
boundary layer (of thickness 0(Re ~ 1 ^ 3 ) where Re is the global Reynolds number for the 
flow). This inner, classical boundary layer develops a Goldstein (1948) singularity at a 
position downstream of the position of maximum channel constriction. However, Smith 
and Daniels demonstrated that this singularity could be removed, without any upstream 
influence being induced, and their analysis was able to describe the reattachment of the 
boundary layer at a distance 0(/i 3 ) downstream. 

Subsequently Hall and Bennett (1986) considered the linear instability of this flow, 
in the general context of the stability of both steady and unsteady interactive boundary 
layer flows, at stream-wise locations before the Goldstein singularity develops. This study 
demonstrated that the flow is unstable to Gortler vortices for channels in which the con- 
striction has a region of concave curvature, in an interval before the position at which 
separation occurs. In this analysis / rx , where f(x) = hF(x ) with h ^ 1 is the constriction 
profile, effectively plays the role of the Gortler number. 

In the case of external boundary layer flows over a wall of variable curvature Hall 
and Lakin (1988), motivated by the earlier weakly nonlinear theory of Gortler vortices by 
Hall (1982b) and the numerical solution of the fully nonlinear Gortler vortex governing 
equations by Hall (1988), developed a self consistent asymptotic theory for short wave- 
length, large amplitude Gortler vortices. In this theory the Gortler number G scales with 
A; 4 (where k 1 is the vortex wave number). Their analysis demonstrated that the mean 
flow adjusts due to the presence of the vortex so as to render all vortex modes neutrally 
stable. Within an 0(1) region of the boundary layer the mean flow is driven by the vortex 
velocity field and as such has no relation to the unperturbed boundary layer flow. Such a 
situation is similar to that postulated by Malkus (1956) for turbulent flows; Malkus argued 
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that the mean part of a turbulent flow would adjust so as to ensure that all ‘modes’ are 
rendered neutrally stable. 

As pointed out by Hall and Lakin (1988), the modification of the mean flow by the 
nonlinear longitudinal vortex motion raises the question as to what effect this modified 
mean flow has on the onset of separation, in separating flows that can support Gortler 
vortex like disturbances.- One important result of the Hall and Lakin calculation is that 
the skin friction for the vortex induced mean flow can be dramatically increased over that 
of the unperturbed boundary layer (although no quantitative results on this point are 
presented this can be inferred from their figure 6, Hall and Lakin (1988) pp 440-441); we 
would then expect that the onset of separation could be considerably delayed due to the 
presence of the vortex motion. In light of the work of Hall and Bennett (1986) alluded to 
above, we choose to concentrate on the Smith and Daniels (1981) problem of the steady 
flow in a symmetrically constricted channel. 

The procedure adopted for the rest of this paper is as follows. In §2 we develop 
the governing equations for nonlinear Gortler vortices in the 0(h~ l ! 2 ) classical boundary 
layer of Smith and Daniels together with a discussion of the relevant results from Smith 
(1982) concerning the onset of separation as well as the those points from the linear theory 
of Hall and Bennett (1986) pertinent to the present investigation. In §3 we extend the 
linear theory of Hall and Bennett (1986) to the fully nonlinear regime for short wavelength 
vortices considered by Hall and Lakin (1988). In §4 we present the numerical solution 
of the vortex induced mean flow equations. Finally in §5 a discussion of our results is 
presented; particular emphasis is given to the effect that the presence of the nonlinear 
vortex structure has upon the onset of separation. 

2. Formulation of the governing equations 

Consider the steady motion of an incompressible fluid of density p and kinematic 
viscosity v in a symmetrically constricted channel. Let L denote the typical length scale of 
the constriction and let Uoo denote the typical free-stream speed sufficiently far upstream 
of the constriction (see figure 1). We define the Reynolds number Re by Re = U^L/v, 
and throughout this work we will assume that the flow is such that the Reynolds number 
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is asymptotically large. The steady three-dimensional Navier-Stokes equations in nondi- 
mensional form are 


du dv dw 

dx + dy + dz ’ 

du du du 

u— + v— + w-z- 
dx dy dz 

dv dv dv 
u— + v— +w — 
dx dy dz 

dw dw dw 

u— + v— +w — 
dx dy dz 


dp _j / d 2 u d 2 u d 2 u 
dx + \ dx 2 dy 2 dz 2 



( 2 . 1 ) 


where ( x , y, z) are the usual Cartesian co-ordinates, nondimensionalized with respect to 
L, (u, v, w ) are the corresponding velocity components, nondimensionalized with respect 
to the upstream velocity Uoc and p is the pressure, nondimensionalized with respect to 
pU^. For definiteness we will assume, following Smith (1982), sufficiently far upstream of 
the channel constriction the flow is Poiseuille 


u 


y-y 




2x 

Re* 


x 


-oo. 


As shown by Smith (1982), with a constriction length L — 0(1), the critical height of the 
constriction is 0(Re~ i ^ 3 ) (Smith’s ‘strong constriction’). We then take the constriction 
y = hRe~ 1 ^ 3 F(x) at the lower wall (a corresponding constriction y = 1 — hRt~ 1 / 3 F(x) 
is also present at the upper wall, however here we are predominantly concerned with the 
boundary layer flow in the vicinity of the lower wall and as such we will restrict our 
attention to the lower boundary); we will subsequently take the limit h 1 (this being 
the appropriate asymptotic regime of the Smith and Daniels (1981) investigation) however 
we initially consider 0 < h < oo. The flow develops a viscous wall layer of thickness 
0(Re~ 1 / 3 ), in which the velocity field and pressure expand as (Smith (1982)) 

(u,v,w,p) = (Re-'l 3 U,Re- 2 l z V,Re-V z W,Re- 2l3 P Q + Re~ 4/3 P 1 ), 
x = 0(1), y = Re~ 1 / 3 Y, z^Re~^ 3 Z, 


( 2 . 2 ) 
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where P 0 is a function of x while P\ is a function of x,Y,Z. The governing equations in 
the viscous wall layer are then given by 

dU dV dW _ 
dx + dY + dZ ~ ’ 

TT dU , r dU dU dP 0 d 2 U d 2 U 

dx ^ dY^ dZ dx dY 2 + dZ 2 ’ 

TT dV dV dV dP\ d 2 V d 2 V 

dx dY dZ dY dY 2 dZ 2 ’ 

TT dW „dW TTT dW dP i d 2 W d 2 W 

U b V b W — n 1 

dx dY dZ dZ dY 2 dZ 2 ’ 

together with the boundary conditions 

(u,v,w,p) — > (y, 0,0,0) as x — > — oo, 

u = v = w = 0 on Y = hF(x), 

U~Y oo. 

Here we have anticipated the well known result that the Gortler vortex wavelength is scaled 
on the boundary layer thickness. 

In the asymptotic limit h 1, Smith and Daniels (1981) have demonstrated that the 
flow described above develops a classical boundary layer of thickness 0(/i -1 / 2 ) attached 
to the hump. In the absence of any vortex component the appropriate scalings are found 
to be 

U = hu, V = h' l ' 1 v, P Q = h 2 p ( = h' /2 Y . 


(2.3) 


From the analysis of Smith and Daniels (1981) the governing equations are given by 

(2.5a, b) 

ux uq 

together with the boundary conditions 


_ du _du d 2 u du dv 

U Ih +V dZ =z x + dt 2 ' lh + dt =0 ' 


u = v = 0, £ = 0, u —* F £ — > oo, 


(2.6a, b ) 


and the upstream conditions 


u,v — > 0 x —> — oo; 


(2.6c) 
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(here we are assuming F(x ) — ► 0 as x — > — oc and we have made use of the Prandtl 
transformation). We note that in (2.5a) the pressure gradient term —FF X arises from the 
outer flow solutions (the reader is referred to Smith and Daniels (1981) for full details). 
Subsequently we will assume that the constriction profile is such that F(x) = 0 for x < 0 
so that the upstream condition (2.6c) can be replaced by 

u,v = 0 x < 0. (2.6<f) 


To derive the governing equations for the Gortler vortex velocity field, which we 
assume is confined to the secondary 0(h 1/2 ) boundary layer, the following scalings are 
employed (see Hall and Bennett (1986)). The x,y and z disturbance velocity components 
scale with the x,y and y basic state velocity components; the y and z variations are on 
the classical boundary layer scale, namely 0(/i -1/2 ), and the pressure perturbations are 
reduced in size until the pressure gradient and viscous terms in the y and z momentum 
equations become comparable. Thus 


£ = h'/ 2 Y, z = h^ 2 Z, 

U = hu, V =h l/2 v, W = h 1/2 w , P 0 = h 2 p, Pi = hp. 


(2.7) 


The governing equations are then given by 


dz 2 ’ 


du dv dw 

Ox o£ Oz 

du du du d 2 u d 2 u 

u di + "af + ” Tz = FF ' + ee + a?’ 

dv dv dv ,, /2 „ 2 dp , d 2 v , d 2 v 

U di + V M + W 9^ + ““ =- af + ae + a?’ 

dw dw dw dp d 2 w d 2 w 

U d^ + V d( +W dl"^ + W + 

together with the boundary conditions 

v = w = 0 £ = 0, 

► P, V£,w —* 0 ( — > oo, 

u, v, w = 0 x < 0, 

where we have made use of the Prandtl transformation. 


u 


u 


(2.8a - d) 


(2.9 a — c) 
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Our concern is with the solution of the nonlinear system (2.8,9) in the limit of short 
spanwise wavelength, djdz 1. Before proceeding with a description of the solution of 
these equations we make some comments concerning the unperturbed system (2.5,6) and 
recall the relevant details from the linear theory of Hall and Bennett (1986). 

In practice the solution of the system (2.5,6) is a numerical task; we briefly describe 
the numerical technique used to solve this system. In order to obtain solutions of (2.5a,c) 
to within the numerical accuracy required for the subsequent calculations of §4 it was 
found to be advantageous to solve a modified system obtained from (2.5) by writing 


u = Fu 0 , v = —F x £ + n o, 


where (uq,uo) are governed by 


F x ul + Fu 0 ^p + (v 0 - F S Q^£- = F x + 


du 0 


du 0 


dx 


dt 


d 2 u o 
d( 2 ’ 


_ _<9u 0 dv o 

F,(uo - 1) + F— + — = 0, 


with the appropriate free-stream boundary condition now given by 


u o 


i 


oo; 


all other boundary conditions remaining unchanged. This modified system was discetized 
according to 


- p(«„ + , - 2fi„ + M + F , (1 - u 2 „) 


(u n fjt(n)(Un+ 1 U n — j) 


^n+1 ^n— 1 F(u n U n ) 

2 h = e 


2 h 

- F x (u n - 1), 


(2.10a, b ) 


where a tilde denotes a quantity evaluated at x + e and the index n refers to a quantity 
evaluated at the grid point £ = nh ; we explicitly impose the boundary conditions (2.6a,b,d). 
Given u n , v n and F the tridiagonal system (2.10a) can be solved for the updated value u n \ 
v n is then determined from (2.10b). The solution was marched downstream from x = 0, 
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with the initial conditions (2.6d), iterating at each step until converged values of u and v 
were obtained. For definiteness we will restrict our attention to the hump profile 


1 0 , x < 0 , 

ae -1 (x 4 /16 — x 3 + 4x 2 ), 0 < x < 4, (2.11) 

ai 2 e -1 ^ 16 , x > 4, 

where a. is an amplitude parameter. Consider now the boundary conditions (2.9b). As 
noted by Smith (1982) the boundary layer is accelerated where F x (x ) > 0 (upstream 
of the point of maximum constriction x = r ma x = 4). Downstream of x max the flow is 
decelerated due to the adverse pressure gradient — FF X ] the skin friction then tends to zero, 
the Goldstein singularity is encountered and the boundary layer flow becomes separated. 
In figure 2 we present a plot of the skin friction r(;r) = u^(x,0) versus x for the profile 
(2.11) and the particular choice a = 1.031. We see that r attains a maximum at x « 2.9 
and then decreases until r = 0 at x = x s ~ 4.53 at which point the Goldstein singularity 
will arise and separation will occur. We note that changing the amplitude a has no effect 
on either the position at which r attains its maximum nor the point at which r = 0. Our 
interest lies in the effect that the presence of a fully nonlinear vortex structure upstream 
of imix may have in delaying (or promoting) the onset of separation. 

The linear theory of Hall and Bennett (1986) demonstrates that the boundary layer 
becomes unstable to linear Gortler vortices at the position x* where 


2F xx u 


du 

ae 


1 , 


(**,*•) 


( 2 . 12 ) 


and are confined to a thin layer about the position £ = where 


d 

dt 



= 0 . 




(2.13) 


As noted by Hall and Bennett, for humps with F > 0 (ie channel constriction) uu^ remains 
positive until the Goldstein singularity develops, at x — x s say, and hence for instability we 
require the constriction to be locally concave somewhere in the interval 0 <■ x < x s . The 
profile (2.11) was chosen with this in mind; it is locally concave in the region 0 < x < 1.69. 


7 


As an aside, we note that, from the results of Hall and Bennett, there is a critical 
amplitude a below which instability cannot occur (see expression (3.15) and the preceding 
discussion in Hall and Bennett (1986)). In figure 3 we present a plot of the neutral 
position x* versus the amplitude a for the profile (2.11); we find the critical amplitude to 
be a = a c k 0.84. For a > a c there is a finite interval Xq < x* < x* in which the flow can 
support linearly unstable Gortler modes and this interval increases with increasing a. 

3. Evolution equations for large amplitude Gortler vortices 

We now proceed to develop an asymptotic solution of the governing equations (2.8a- 
d) valid in the short wavelength limit. We assume that the solution (u(r, £), v(x , £)) of the 
classical boundary layer equations (2.5a,b) becomes linearly unstable to Gortler vortices 
at x = x* where x* is to be determined from (2.12). 

In the 0(h~*l 2 ) boundary layer the effective Gortler number G becomes 0(h 3 / 2 ) so 
that, on the basis of the nonlinear theory of Hall and Lakin (1988), we anticipate vortices 
with wavenumber k = 0(h 3 / 8 ) (the scaling G = 0(k 4 ) is that which is appropriate to the 
right-hand-branch of the linear neutral curve, see Hall (1982)). At first sight these scalings 
appear to be in disagreement with those presented in Hall and Bennett (1986); in fact, 
there is no discrepancy when one considers the boundary layer scalings (2.7). The analysis 
of Hall and Lakin demonstrates that in the presence of the nonlinear vortex state the 
flow regime breaks up into three regions; a central core region, in which the mean flow is 
induced by the vortex velocity field, which is bounded above and below by thin shear layers 
in which the vortex amplitude is reduced to zero. Outside the region of vortex activity the 
mean flow is governed by the unperturbed boundary layer equations (see figure 4). We now 
proceed to develop the corresponding structure for the flow described in §2; we emphasize 
that the nonlinear vortex state will be confined wholly within the 0{h ~ ] / 2 ) boundary 
layer. Hence, we are anticipating that the only effect on the onset of separation due to the 
presence of the nonlinear vortex state is that the position at which separation occurs will 
be changed; the asymptotic analysis of Smith and Daniels concerning the removal of the 
Goldstein singularity still remains valid. As the analysis mirrors that of Hall and Lakin we 
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will give only the salient features of that calculation (the reader is referred to that paper 
for full details). 

To proceed we first consider the core region, in which the vortex amplitude is chosen 
so as to have an 0(1) effect on the mean flow in this region. The appropriate expansions 
in the core region are then given by 

u = uq -f- h 3 l t/oi cos h 3 ^ s z + • • • 1 , 

r ^ M 

v = v 0 + h 3 ^ s I Voi cos h 3 / s z + ■ , 

together with similar expansions for w and p; here • • • denotes terms which do not enter into 
the subsequent analysis. Substituting the expansions (3.1) into the governing equations 
gives, for the vortex components 

Uqi + uq^Vqi = 0, Vq\ + 2uqF xx Uo\ = 0, (3.2a, 6) 


which for consistency requires 

2F xx u 0 u h — 1 . 

Integrating (3.3) yields, for the vortex induced mean flow 


u 0 — 


A(x) + £ 


(3.3) 


(3.4a) 


while from the continuity equation we have 


vo 


— a (* ) -W^ + ¥(^r 1 ) 


3/2 


(3.46) 


here A(x) and B(x) are to be determined. Thus, in the core region the mean flow is found 
from a solvability condition on the vortex amplitude equations. At this stage the vortex 
amplitude Voi remains undetermined; this is found from the mean flow component of the 
first momentum equation. Thus 


_ du 0 _ du 0 d 2 uo -- 

U ° Ox + u ° d£ di 2 x ~ 2 di V d( V01 1 ’ 


1 d ( du 0 y2 
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which making use of (3.4) can be integrated to give 


C(x) + 


4 y/F xx (A + <!;) 


Co 2 , = -B, 


r A + ( F„,(A + (f 


6 *?, 


2-y^ F XX (A + £) 


-FF X (. (3.5) 


If we let £j(x) (j = 1,2) denote the outer limits of the region of vortex activity we then 
have from (3.5), upon eliminating C(x), a relation between and £ 2 , namely 

A + £1 _ F x XX (A + £1 ) 2 1 

Fxx 6 F% x 2y/F^JA 4 - £ 1 ) 

_ FxxxjA^ 4~ ^ 2 ) 2 1 

2\^~Fxx(A + £ 2 ) 



FF x i , = -B 
FF X &. 


A + £2 


(3.6) 


Thus, the vortex amplitude in the core region is reduced to zero as £ — ► £i > 2 ; however we 
require thin shear layers centred on £ 1,2 in order to smooth out the algebraically decaying 
vortex component. We note that at this point of the analysis we are unable to determine 
the functions A(x ), B(x :), £i(x) and £ 2 ( 3 :)* 

We now consider the solution in the shear layers. As noted by Hall and Lakin, the 
thickness of these layers is found by balancing diffusion across the layer with convection 
in the stream-wise direction. In the case under consideration such a balance shows the 
thickness of these layers to be 0(h ~ 1//4 ). We then write 


n = fc 1/4 (e - 6), 


and in the shear layer we replace djdx and d/d£ by d/dx — h l / 4 £ f 2 d/dr} and h l ^ 4 djdi^ 
respectively. We will subsequently restrict our attention to the shear layer centred on £ 2 ; 
an identical analysis holds in the layer centered on £ = £ 1 . 

From the core solutions (3.4-5) the appropriate expansions in this shear layer are 
found to be 

u = Uq + /z" 1/4 ui + /i _1 / 2 u 2 + * - * + h~ A t 2 |(t/oi + h~ l t 4 u\\ + • • •) cos /i 3 / 8 z + • • • J , 
v = Vo + h-^% 4 f |(«oi + h~^ 4 v u 4 )cos h z ^z 4 j , 
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together with similar expansions for w and p. Substitution of these expansions into the 
governing equations we find that, after some trivial manipulations, the leading order mean 
flow terms are given by 


u o 


l~A + £2 = 

— = , ui = 


XX 


vq = -B{x) - Aj 


2\/ F XX {A + £ 2 ) 


I A + £2 -Fxxx 

m I r* 


(3.7a, c) 


where we have made use of matching with the core solutions in the limit 77 — * — 00 . The 
governing equations for the leading order vortex amplitudes (t/o 1 , ?! o 1 ) are found to be 
consistent whereas the equations for (uu,un) are consistent provided that 

d 2 V(n S(x)pv 01 1 3 , . 

— = 5 


which gives the governing equation for the vortex amplitude t>oi • Here 

S = -2B- { ^ 1 } 3/2 + -tt T - + 

3 { F xx ) A + & 

and fi is a funcion of x , the precise nature of which is not required here. This equation is 
a modified form of the second Painleve transcendent and has been shown to have a unique 
solutiSn, Hastings and McLeod (1980), with the asymptotic behaviour 


Vqi ~ Sr 7, p — > -00, u 0 i 0 , p -* 00. 

It then follows that the vortex is constrained to lie below £ = £2 ( a similar analysis of the 
shear layer about £1 shows that the vortex is constrained to lie above f = £ 1 ) so that the 
nonlinear vortex is confined to the 0(1) core region. We note from (3.7a— c) that the mean 
flow in the shear layer is essentially unaltered by the presence of the vortex (in fact, the 
first two terms in the stream-wise mean flow expansion in this layer are obtained from 
the Taylor series expansion of the core flow about (2 )• Hence the mean flow outside the 
core region must have u, and v, in the limits £ — ► defined by the core solution 

evaluated at £ = £i j2 respectively. 
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Outside the region of vortex activity there is only a mean velocity field and we write 


u — u + 0(h 1 / 2 ), v = v + 0(h J / 2 ), 


where (u,v), defined in (0,£i) and (£ 2 , 00 ), satisfy 


_du _du d 2 u __ 

U fc +V d( = dp + x ' 


du dv 

di + oi = °- 


This system is to be solved subject to 


u = u = 0, £ = 0, u — » F, £ —y 00 , 


(3.8) 


(3.9a, b) 


and 


v(£j) - ~B(x) - A 2 


U(($j) = 


2 y/F^A + b)’ 


A ~h F xxx 

F xx 3 


(^r) 


3/2 


0 = 1 , 2 ) 


(3.9c -e) 


The system (3.8,9) together with (3.6) then specifies a free-boundary problem for £i(x), 
£ 2(^0 and the functions A(x), B(x) together with the skin friction u^(x,0). The solution 
of this system is a full numerical problem (we note that the similarity reduction available 
to Hall and Lakin is not applicable here). Before proceeding with a discussion of the 
marching procedure used to solve this system we note that, in order to solve (3.8,9), it is 
necessary to specify an ‘initial state 7 for the mean velocity field from which the solution 
can be obtained in a systematic way by marching downstream. This requires finding an 
asymptotic form for the solution close to the position x* at which the original boundary 
layer becomes unstable. 

We assume that the unperturbed boundary layer flow first becomes unstable at the 
position x = x* and that the instability originates in a thin layer centred on £ = £*. We 
define X = x — x* 1 together with the similarity variable 77 = (£ — £*)/X 1 / 2 and expand 
£1 and £2 in the form 


6 =r-* i/2 e + op 0 , 
6 = C+x 1 ' 2 ( + 0 (x), 


(3.10) 
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(these expansions follow from the weakly nonlinear theory of Hall (1982b)). Consider first 
the core region; following Hall and Lakin (1988) we express (uo^o) i n the form 

u 0 = u + X 3 / 2 U M (r)) + ---, v 0 = v + XV M (r}) + ---, (3.11a, b) 


where (u,v) are the boundary layer velocity components in the absence of the vortex 
motion. In a small neighbourhood of the neutral location (#*,£*) we have the Taylor 
series expansions for (u,v) and F xx given by 

u = Uoo + X'/ 2 7jUio + X(rj 2 U 20 + «oi) + X 3 ^ 2 rjun + • • • , 
v = woo + X l ^ 2 rjV\Q + X(r] 2 V 2 o + v oi) + X 3 ! 2 r)V\\ + • ■ ■ , (3.12a — c) 

F xx = Fq + XF\ + • • • . 


Substituting (3.11a,b) together with (3.12a-c) into (3.3) and equating like powers of X 
gives 


2Fouooaio = 1, Wjq + 2uqo u 20 — Oj 


where 


Ao 


Um = A 0 r/ - Air/ 3 , Vm — 
aio^oi Ti^io 




+ K, 


(3.13a -c) 


-an - 


Ai = 


U10U20 


uoo Fq Moo 

and V£j is a constant of integration (required to ensure continuity of Vm across the shear 
layers). We note that (3.13a) is automatically satisfied as x* is the neutral position, while 
(3.13b) is satisfied since £* is chosen so that (2.13) holds. 

Following Hall and Lakin (1988) we can show that the vortex amplitude Voi vanishes 
at r] — rj± where 


\ri±f = 


Cq + 2Aq 


(3.14) 


6Ai + uooAo 

with Co an unknown constant to be determined. It then follows from (3.10) and (3.14) 
that 

v- = -L r i+ = c 


To determine r}± it remains to determine the constant Cq. 
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Outside the region of vortex activity we now write 


u 


= u + X 3 / 2 Um(t]) + • • • , v = v + XVm{v) + 


We note that the functions U m and Vm are, respectively odd and even functions of r/ and 
so we will restrict our attention to the region ij > r) + . Here we find that U M ,V M satisfy* 


d 2 UM uoo ( d \ 
di? + 2 \ V dj] 

dV M 3 r/ dU M _ 

+ 2 Um ~ 2~W = 


which are to be solved subject to the matching conditions 

Um = A 0 r/+ - Air/® , U' M = \ 0 - 3\iT]+, 

and 

Um — > 0 rj — ► oo. 

The appropriate solutions are found to be* 

JJ = (W - Xivl)exp {-u 00 (? 7 2 - r?j.)/8} 17(7/2, y^r/) 
M 17(7/2, v^sr, 17+) 




Vm = ~ ^ / Z7 M dr;, 


7 

1 


(3.15) 


(3.16a) 


(3.166) 


(3.18a, 6 ) 


where U(a,x ) is the parabolic cylinder function. From (3.17a) and (3.18a) we then have 
the eigen-relation* for 77 + 


_ / ^oo f (9/2, 

' T\ U{7/2 ty /*?r, + ) 


+ 


«oo I A 0 - 3A it ]1 
*1+ } = 


(3.19) 


Ao»?+ - Ai^ ’ 

which then determines the constant Co in (3.14). 

With the above asymptotic solution, valid near x *, we can construct the composite 
velocity field together with initial values for the function A(x), B(x), £1 (z) and £ 2 (®) which 
will then be used to start the marching procedure to be described in the next section. 


t 


t: 


Note the typographical error in the corresponding expressions of Hall and Lakin (1988). 


bid. 
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4. Numerical solution of the free boundary problem 

The numerical scheme used to solve the free boundary problem (3.8-9) is based on that 
used by Hall and Lakin (1988); for this reason we only give the relevant details here. We 
assume that a solution is known for x* < x < x (which in practice must be constructed 
from the asymptotic solution for the initial development of the nonlinear vortex state 
presented above). The scheme then advances this ‘initial’ solution downstream to x + e for 
some sufficiently small value of e. 

In the course of our calculations it was found appropriate to solve a modified system 
of the form outlined in §2. Thus we define 


u = Fu , v = v — F x £. 

The numerical scheme then proceeds by defining a new variable i/’ by 

r 


^ = 




j = l, o<£<£„ j = 2 , 6<£<°o. 


The modified form of the system (3.12) is now written as 


du 1 d 2 


Fu— - 


u 


dx £ 2 cb/’ 2 


= F x (l-u 2 )- 


(v - F x tj<f>) du Fi) du 


+ 


X \)U‘ 


0 d 6 r dr 


du t f ^ M ^ ( ipFQ du 

u) F dx + d^J' 


(4.1a, b ) 


where j = 1 for 0 < 0 < 1 and j — 2 for 1 < -0 < oo. This system is to be solved on (0, oo) 
in terms of ij) with boundary conditions 


u = u = 0, = 0, u — ► 1 oo, 


(4-2) 


together with the ‘discontinuous’ conditions at xp — 1 


1 A + tj 


U± F \ F x 
v± = —B — A x 


u v ± 


6 


I A i \j Fxxx 

T n 


2F v 'F rx (^ + e i )’ 


(4.3a, b) 
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Here the ± signs corresponds to the limits -0 — > l -1 ", -0 — » 1 — while the index j = 1 is 
associated with the c — ’ sign and j = 2 with the ‘ + ’ sign. The system is then closed by the 
condition (3.6). 

The solution of (4.1a) is advanced for ip in the range 0 < ip < 1 by using the scheme 


Fll n U n ^2^2 (^ n +l 2 U n U n — i ) — Fu n (F £ ( 1 u n ) 

2£~h Fx^jlpn Flpn^jUfi^ (u n -j-] U n — j), 


(4.4) 


where h is the vertical grid spacing, a tilde denotes that the quantity is evaluated at 
x + e while an index n indicates a quantity evaluated at the grid point xj) = nh. In order 
to solve the tridiagonal system (4.4) we must make an initial guess for £1 and then set 
£i = (£1 £i)/ 6 - solving (4.4) we ensure that the boundary condition on u at xp = 0 

and (4.3b), with j = 1, are satisfied. The solution of (4.4) then yields u n . We discretize 
the modified continuity equation (4.1b) as 


^n+l 1 

2 h 


— 




+ 


Flpntj (Un+ 1 ~ U n - 1 ) \ 

2 h j 


(4.5) 


which, given u n from (4.4), allows us to determine v n at the position x + e for 0 < ip < 1. 
As the equation for v is first order in xp only the boundary condition on v at ip = 0 can 
be satisfied. Thus, although the solution of (4.1), for ip in the range 0 < ip < 1, has been 
obtained at x + t the additional boundary conditions (4.3a,c), with j = 1, have not been 
satisfied. These conditions can now be used to obtain a value for A and an improved value 
for by writing them in the form 


€i =F £ F xt ul-A, 


v- + B — 


{A — A) A + 


F 

L XI 


+ 


A + 


k ) 


3/2 


+ F x £ i, 


(4.6a, b ) 


Here B is the current guess of B at x + e. This scheme can be applied in a similar manner 
to the region ip > 1. The equation (4.4) is now solved for u subject to u — > 1 as ip — » oo 


16 



together with (4.3b) for j = 2, and equation (4.5) is solved for v subject to (4.3c) with 
j = 2. Writing (4.3a) and (3.6) in the form 



v 


6 = F 2 F xx u\ - A, 


^4 + £2 


F. 


6 J?, 


{(i + fi) 2 -(i + 6) 2 } 


+ — b = I , 1 - , 1 .■ > + fM> - 6), 

2y/F xx [ yj a + ii yj A + I 2 ) 

(4.6c, d) 

the first of which then determines a new value for £2 while the second equation yields a 
value for B. We then have values of u, v, A , B, £1 and ( 2^1 + e. This procedure is then 
continued, making use equations (4.6a-d), using Newton iteration on A,B,£\ and £2 until 
converged values of these quantities are obtained. 

For definiteness, our numerical calculations were restricted to the profile (2.11). The 
first stage in the numerical solution of the free boundary problem (3.8,9) involves deter- 
mining the neutral position from (2.12,13). This was achieved using the numerical scheme 
outlined in §2; with the initial conditions « = v = 0atx=0 the solution was marched 
downstream, at each stream— wise location was determined from (2.13), until the point of 

neutral stability x* , determined from (2.12), was reached. At this point an initial guess for 
^(x), £ 2 ( 20 , A{x) and B(i), for x = x* + 0.001, was obtained using the procedure outlined 
at the end of §3. The initial condition required for the velocity field at x, namely u ( .? , £ ) 
and u(x, £), were obtained using the scheme (2.10) and these values were then interpolated, 
using cubic spline interpolation in order to ensure continuous second derivatives, to yield 
u and v on a uniform grid in ip i n the intervals [0, 1) and (l,t/>oo). A similar interpolation 
procedure was used on the vortex induced mean flow corrections U m ,V m thus allowing the 
composite velocity field to be computed. The composite velocity field so obtained was then 
used as the initial conditions for the marching scheme described above. The solution was 
marched downstream until such point as the limit of the region of vortex activity tended 
to zero (ie at the point xp where £1 (xp) = £ 2 ( 2 ^), see figure 4). At this point the velocity 
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field was then used as the initial condition for the scheme (2.10) which was then used to 
march downstream until such point at which reversed flow is first encountered. 

5, Conclusions 

Before proceeding with a discussion of our results we make some comments concerning 
the stability characteristics of our numerical scheme. Initially the normal step-size h was 
fixed at h = 0.01 with xftoo — 16 thus yielding 100 and 1500 normal grid points in the 
intervals 0 < rfi < 1 and 1 < < ^oo, respectively. With h fixed, various calculations 

were performed with different values of the streamwise step-length e; from e = 0.001 to 
e = 0.00025. For the range of e considered the position x/r, at which the vortex region 
vanishes, was found to vary by, typically, 0.9% over the value obtained for the smallest 
value of e considered. Similarly, the normal position was found to vary by 0.7% over 
the range of e considered. We deem our numerical scheme to be stable and for e ~ 0.00025 
suitable accuracy has been attained. All of the subsequent results were obtained for the 
choice of parameters h = 0.01, e = 0.00025 and — 16. 

Consider figure 5, in which we present a plot of the positions (x) and &( x ) bounding 
the region of vortex activity for the constriction amplitudes a = 1.0304 and a = 1.0923. 
As anticipated from the results of Hall and Lakin (1988) we see that the region of vortex 
activity initially thickens until such a point at which the maximum normal extent of this 
region is achieved. As x is further increased, the region of vortex activity constricts until 
the point x? at which £\(xp) = £ 2 { x f)- For x > X F the flow can no longer support the 
large amplitude vortex state described above. 

The effect of increasing the amplitude a of the constriction is evident from a com- 
parison of figures 5a, b. The minimum values of £1 in figures 5a, b are 2.19 x 10~ 2 and 
1.08 x 10 2 respectively whereas the maximum values of £2 are 1.384 and 1.752 respec- 
tively. Thus the effect of increasing a is twofold; firstly the streamwise extent of the vortex 
region is increased with increasing a (a consequence of the change in the neutral positions 
as represented in figure 3) and secondly, the normal extent of the vortex region increases 
with increasing a, through both the increase in the maximum value of £2 and a decrease 
in the minimum value of £ 1 . 
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Turning now to the question as to what effect the nonlinear vortex state, present 
in the interval (x*,X/r), has upon the onset of separation we present, in figures 6a-b, 
plots of the skin friction in the presence of the vortex (broken line) and, for comparison, 
the corresponding skin friction for the unperturbed boundary layer (solid line), for the 
representative cases a = 1.0304 and a = 1.0923, Here the skin friction is defined according 
to 

nf^ du ( X ’°) _ F ( X ) 9u(x,0) fRU 

x \ x ) > / \ o / ? \*/ 

d( 6W dip 

where ip is the modified independent variable and u is the modified dependent variable used 
in the numerical scheme described above. From a computational standpoint we replace 


(5.1) by 


F(x ) du(x,0) 

600 dip 


{iu(x, h) - u(x, 2 h)} + 0(h 2 ). 

2 hc,\(x) 


(52) 


The effect on the boundary layer due to the vortex is readily seen through the dramatic 
increase in the skin friction (over that for the unperturbed state). In the presence of 
the vortex the skin friction now develops two maxima, the first being due to the vortex, 
whereas the second is seen to coincide with that for the unperturbed state. The ratio of 
these maxima is seen to increase with increasing a; for the case of figure 6a this is found 
to be 73% while for figure 6b is 93%. However, in both figures 6a, b we see that the vortex 
induced skin friction subsequently decreases to below the value found for the unperturbed 
state followed by an increase, due to the presence of the favourable pressure gradient, until 
such point as the second maxima is attained (at which point the difference between the 
vortex induced state and the unperturbed state is barely discernible). The skin friction 
now decreases until the point of separation is encountered, the location of which is identical 
for both the unperturbed state and the case in which the flow supports a nonlinear vortex. 

We are then given to conclude that the presence of a nonlinear vortex, in a region 
of concave curvature upstream of the position of maximum channel constriction, has no 
discernible effect upon the position at which separation is encountered. This result is 
somewhat surprising in light of the comments of Hall and Lakin (1988) (pp 443) that the 
decay of the vortices does not allow the boundary layer to return to its undisturbed state. 
In the present case it would appear, by the position at which separation is reached, that 
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the boundary layer does, in fact, return to its unperturbed state. In figure 7 we present 
a plot of the streamwise velocity, for the case of the unperturbed boundary layer and the 
case when the vortices are present, at the position xf- 

In the analysis presented we have ignored two other possible mechanisms for linear in- 
stability and subsequent nonlinear growth, namely Tollmien-Schlichting waves and inviscid 
Rayleigh modes (see Tutty and Cowley (1986) for a discussion of the Rayleigh instability 
in interactive boundary layer flows). The nonlinear development of these modes, together 
with the possible importance of the vortex-wave interaction mechanisms described by Hall 
and Smith (1988), could have a dramatic effect on the results and conclusions of this work. 

Finally, for the (somewhat more realistic) case of moderately constricted channels, for 
which h = 0(1), the effective Gortler number is now an 0(1) quantity and we would expect, 
on the basis of the linear theory of Hall (1983) for vortices with an 0(1) wavenumber, that 
the flow will support linearly unstable Gortler modes. The extension into the nonlinear 
regime then involves the numerical solution of the full nonlinear governing equations as 
was carried out by Hall (1988) for the external boundary layer problem. In the present 
case the governing equations show distinct similarities with those derived by Denier and 
Hall (1992) to describe the nonlinear evolution of the most unstable Gortler vortex mode 
the solution of which resulted in a separated flow. In the context of an 0(1) constricted 
channel, the resolution of this problem would require the numerical solution of the full 
nonlinear governing equations. This is currently the subject of ongoing research. 
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upstream, flow 



Figure 1 . A sketch of the flow configuration under consideration. 
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Figure 2. Plot of skin friction r = ^0,0) versus x for the profile (3.3). 
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Figure 6b. See figure 6a for legend. 
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Figure 7. Streamwise velocity profiles for x = xp. The vortex induced mean flow (broken 
line) and the unperturbed mean flow (solid line). 
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